Spin configurations of carbon nanotube in a nonuniform external potential 
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We study, theoretically, the ground state spin of a carbon nanotube in the presence of an external 
potential. We find that when the external potential is applied to a part of the nanotube, its variation 
changes the single electron spectrum significantly. This, in combination with Coulomb repulsion 
and the symmetry properties of a finite length armchair nanotube induces spin flips in the ground 
state when the external potential is varied. We discuss the possible application of our theory to 
recent measurements of Coulomb blocked peaks and their dependence on a weak magnetic field in 
armchair carbon nanotubes. 

PACS numbers 61.48. +c, 73.23.-b, 71.10.-w, 71.24.+q, 75.10.Lp 



Carbon nanotubes have attracted continuous attention 
since their discovery §. Their unique mechanical and 
electrical properties enable to design and build several 
systems, some may be used for practical applications. 
When they are rolled up in a certain way the electron 
states near the Fermi level have a one-dimensional (ID) 
character [9 . The Coulomb interaction between the elec- 
trons in low-dimensional systems, especially in ID, affects 
significantly the physical properties of the electrons. In- 
deed, power law behavior of the tunneling conductance 
was found experimentally j|. Recently, Coulomb block- 
ade peaks were investigated in a single-electron transport 
experiments through short (~ 0.2 — 3 /im) isolated metal- 
lic nanotubes || 0j. 

In the experiments of Tans et al. (I) and Cobden 
et al. H (II) the spin direction of the additional elec- 
trons entering into the tube, as an external gate potential 
is varied, was determined various methods. Within the 
framework of the "orthodox theory" , where effects of the 
Coulomb interaction are modeled as a constant charging 
energy, one expects that the electrons will occupy states 
in pairs of spin up and spin down. The last level is singly 
or doubly occupied depending on the parity of the total 
number of electrons in the system. In that case, succes- 
sive electrons entering the nanotube at a low magnetic 
field should have opposite spins. This behavior was ob- 
served in II, but in I eight successive electrons entered 
the tube with parallel spins |l(|. Shorter parallel spin 
(PS) sequences were reported in Ref. ||. 

The purpose of the present note is to explore one pos- 
sible mechanism for PS sequences to occur. The essential 
feature is the existence of two (or more) families of single- 
electron states in the nanotube, which respond differently 
to the applied gate voltage V g and can cross (or nearly 
cross) at certain values of V g . Electron-electron interac- 
tions can then cause spin flips (i.e., internal transitions) 
which can lead to PS sequences under proper conditions, 
as was noted in I. We point out that this can occur in 
ideal defect-free armchair nanotubes, as well as in nan- 
otubes with sufficient disorder that states at the Fermi 
energy are localized in different regions of the tube. 



For an infinite armchair nanotube, for a given wavevec- 
tor parallel to the tube, there are two energy bands of 
different symmetry, which we designate as B and A, dis- 
tinguishing states which are bonding or anti-bonding for 
adjacent carbon atoms around the circumference of the 
tube. As we show below, levels from the B and A bands 
respond differently to a gate potential which is nonuni- 
form along the tube axis. We assume here ideal reflecting 
boundary conditions at the ends of the nanotube, which 
preserve the symmetry, so that the discrete single particle 
levels from the two bands can actually cross each other as 
the gate potential is varied. (See Fig. [I].) If there is weak 
mixing between bands on reflection from a tube end, or 
due to interactions with the substrate, near-crossings can 
still occur, and our conclusions will be unaltered. 

Assume now, for example, that there are 2N particles 
in the nanotube and that the last B-level below the Fermi 
energy is occupied by two electrons with total spin 0. As 
the nonuniform potential (induced by the gate) is de- 
creased an A-level crosses the highest (doubly) occupied 
-B-level at the point Fl in Fig. |l|a. In analogy to the first 
Hund's rule for partially filled (degenerate) shell in an 
atom, the two electrons form a many-body ground state 
with a maximum possible spin. They occupy the A- and 
the i?-levels near Fl point with total spin one as long as 
the energy difference between the levels is smaller than 
the energy gain due to the electron-electron interaction. 
There is a gain in the interaction energy because: (i) the 
direct interaction between the two electrons in the same 
state is stronger then the direct interaction between the 
electrons in different states; and (ii) electrons occupying 
different levels can minimize their interaction energy by 
maximizing the total spin, thereby allowing a maximally 
antisymmetric coordinate wave function, and gaining the 
negative exchange energy. 

Now, assume that the nanotube initially has 2N — 1 
particles where only the last state is singly occupied so 
that the total spin is 1/2. The 2Nth electron enters the 
system at such a gate potential that the ground state 
with total spin one is formed. Then, as the nonuniform 
potential is changed further a spin flip occurs and the 
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spin of the system drops to zero before the (27V + l)th 
electron enters into the nanotube. In this scenario both 
the 27Vth and the (27V+ l)th electrons raise the total spin 
of the system by 1/2. But due to the internal spin flip the 
total spin of the system is not larger then one. The entry 
points of such a scenario are depicted in Fig. [l)b by black 
triangles. We show below that a necessary condition for 
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FIG. 1. (a.) (Top) Evolution of the discrete sin- 
gle-particle levels of an armchair nanotube when a nonuniform 
electrostatic potential <j>{x) on the tube is varied (by chang- 
ing the voltage, V g , on the gate electrode). The dashed lines 
represent the (anti-bonding) A-band, while the dashed-dot 
lines represent levels in the (bonding) B-band. The descend- 
ing solid-line represents the highest occupied level for a tube 
with a fixed even number 27V electrons. Heavy portions of 
the curve show where the highest state is doubly occupied 
(S = 0), while lighter regions near the crossing points Fl, F2, 
F3, show where two states are singly occupied (S = 1) due 
to the electron-electron interaction. Similar S = 1 states will 
occur for 27V + 2 electrons near the crossing points fl and 
f2. For odd-electron number, however, such as 27V + 3, no 
spin-flips occur and the highest filled states follows the tra- 
jectory indicated by the heavy dashed and dot-dashed lines. 
Note that S = 1 states do not occur near the crossing points 
Nl and N2. (b.) (Bottom) The spin S of the system for 
fixed numbers of the electrons in the system when V s is var- 
ied. Electrons added at the black triangles would have spins 
+1/2 while those added at the white triangles would have 
alternating ±1/2 spins. 

the change of the relative positions of the B- and A- levels 
is the presence of a nonuniform gate potential along the 
nanotube. In II Chromium-Gold contacts are deposited 
on the top of a rope with nanotubes, and probably break 
it. As a result, the only active part of the nanotube 
is between the electrodes and is exposed to an almost 
uniform gate potential. (Also, there is likely to be strong 



mixing between the A and B bands at the damaged tube- 
ends.) By contrast, in I, the nanotube lies on the top 
of a Si/SiC>2 substrate with two thick Al/Pt electrodes, 
and it was suggested || that the nanotube is unbroken 
in that area. The length of the conducting part of the 
nanotube (L ~ 3 fim) is then larger than the gap between 
the electrodes (d ~ 0.2 /im). The contact electrodes 
screen the potential applied to the gate, so that only the 
nanotube section above the electrode gap is affected by 
the gate. Thus, in I, the external potential due to the 
gate is not uniform. Relative motion of the A and B 
bands is thus possible, which can lead to spin flips. 

We turn now to a description of our model calculation. 
The energy due to the interaction between the electrons 
is calculated within the Hartree-Fock theory. We assume 
that the electrons occupy the single particle states in 
the presence of a self-consistent external potential on the 
tube <f>(x), where x is the coordinate along the nanotube. 
The total ground state energy in this approximation is 
given by 
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where a discrete level e/ Al [0] corresponds to the Zth 
eigenenergy in the band fi — A, B of the nanotube in the 
presence of the space-dependent potential. (The band 
energies are defined such that the Fermi level is zero 
for a neutral nanotube with = 0.) The occupation 
of the state \lfxcr}, is ni^a — or 1 and 7V Mcr = J2i n l/j.aj 
where a =\,\. The term proportional to magnetic field 
h = g/iBH/2 is the Zeeman energy, and Ni on is the num- 
ber of carbon atoms. The constant Nq is a continuous 
parameter, dependent on the work functions of the nan- 
otube and the contact electrodes, which we choose so that 
the chemical potential is zero in the equilibrium state 
when the contacts and gate electrodes are grounded. The 
values of the charging energy U c , the excess interaction 
of two electrons occupying the same level 6U, and the 
exchange parameter J are found by considering the real 
electron-electron interaction in the nanotube. 

Generically, the interaction contains a long-range part, 
which influences mainly U c , and a short-range part due 
to the local interaction when two electrons occupy a p- 
orbital of the same carbon atom. The amplitude of this 
Hubbard- like term is approximately Uh ~ 15 cV [ pM| . 
The behavior of the long-range part is sensitive to the 
screening by the metallic electrodes. The interaction in 
the section of the nanotube, which is on the metallic elec- 
trodes is screened by the image charges in the metal and 
falls off like R 2 /r 3 for distances r much larger then the 
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nanotube radius R. We find the following values for the 
interaction parameters Jl2| : 

3.94 [5.04] > U c /A> 1.14 [1.34], 
SU/A = 0.11 [0.22], J/ A > 0.22 [0.44], (2) 

for (10, 10) [(5, 5)] armchair nanotubes, respectively, 
where A = itTwf/L with vf = 8.1 x 10 5 m/s is the level 
spacing at the Fermi energy in each band. The lower 
bound on U c is for the screened interaction and the up- 
per bound is for the unscreened interaction. (Higher or- 
der correlation effects between the electrons may renor- 
malize the parameters of the problem [[TJJ.) We have 
checked numerically that the long-range part of the in- 
teraction does not influence J and SU significantly and 
that the effect of a nonuniform potential, which changes 
the self-consistent wave functions, is also negligible for 
the parameter regions we consider below. 

In order to find the evolution of the levels near the 
Fermi energy when the potential on the gate electrode 
V g is changed, we denote the self-consistent nonuniform 
potential on the nanotube by 4>( x )> an( i the energy of the 
last occupied level in band fj, by e M . In the Thomas-Fermi 
approximation the total number of states in band /x up 
to energy e is: 

M^[e,4>]= f de' f dx Pll [e' + <f>(x)] (3) 

J — oo J 

where p M (e) is the density of states in band /i. An applica- 
tion of a potential V g on the gate electrode changes both 
4>{x) and e M but does not change the number A/^ej^, cf>) 
which by definition equals I. Thus for a fixed I we find 
(by formally expanding Af^ with respect to both its ar- 
guments) 

de ifJl _ / dxp^e^ + (j> (x)]d(j)(x) / dV z 
dV g ~ fdxp^ + M*)} 

where 4>o(x) is the self-consistent potential on the tube 
for V g = 0. To model the situation in I we take <f>o{x) 
to be a constant <fio — 0.3V in the regions on top of 
the contacts but zero over the gap, while d(j)(x)/dV g 
is equal to a constant c over the gap and zero in the 
regions over the contacts. For the tight-binding-band 
approximation to an armchair nanotube we may write: 

Pfl (e) cx l/^/l-[(e-^)/(2 7 )] 2 where e° - -e° = 7 is 
the position of /z-band center, and 7 = 2.7 eV is a quar- 
ter of the width of the (extended) bands in the nanotube. 
Thus for <^o^7 we find: 
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where +[— ] refers to the A[B] band. The value of the 
constant c is mainly determined by the ratio of d to the 
distance of the nanotube from the gate electrode. Based 



on a crude model of the geometry in I, we estimate that 
for V g — ±AV the potential on the tube is of order ±0.1 
which gives c « 0.02 JTJ]. 

Using the picture described above, we find that £(^[1^] 
in Eq. ([!]) is given up to additive constants by 

e la [4>] =lA-cVg{a + [3), e lb [4>] = I A - cV e [3, (6) 

where a rj 0.005 and [3 w 0.065 for d/L = 1/15. The 
slope of the levels from A-band (dashed lines in Fig. |l]a) 
is a + (3 and the slope of the levels from S-band (dashed- 
doted lines) is (3. 

Having crude estimations for the parameters of the 
model (|lj) we can find a set of conditions under which 
a sequence of consecutive electrons will enter into the 
nanotube with PS. Fig. [j]b shows the evolution of the 
ground state spin with a fixed number of particles when 
the gate potential increases. The width AV where the 
ground state has spin one depends on the values of a and 
(3 as well as on the interaction parameters 5U and J. 

In order to find an example of sufficient conditions to 
have a PS sequence when electrons enter into the nan- 
otube we first set J = and SU <§; A. For a fixed even 
number 2N of particles we find, after analyzing the evo- 
lution of the highest occupied level (bold line in Fig. pi) , 
that a spin flip occurs whenever an empty A-lcvel crosses 
an occupied i?-level. These crossing points are marked 
by Fl, F2, F3. Notice that if point Fl occurs at a value 
cVg = 0, then fl and 7V1 occur at cV g = A/a, point 
F2 occurs at cV g = 2 A/a, etc. We find that the follow- 
ing conditions are sufficient for a PS sequence: (i) Nq 
is such that the initial Fermi energy is near the crossing 
point Fl; (ii) the (2iV + 2)th electron enters at SV g corre- 
sponding to points fl, /2, ... at which spin flips occur in 
the case of 27V + 2 electrons in the system; and (iii) the 
(2iV + l)th electron does not enters at cV g corresponding 
to the points F2,F3, .... We note that it is possible to 
find other sufficient conditions for a PS sequence, here 
we give one example. Since, with the assumptions J = 
and 5U -€i A, electron entries occur when the distance 
of the line Fl — ► fl — * F2 — ► /2 . . . from the Fermi level 
is equal to the charging energy U c , we find that (i), (ii) 
and (iii) are fulfilled when 

(ii) = (N a - k) mod 4 = 0, 

(ii) U*/A = (1/2) [n + W/a)(2n + 1)] , 

(iii) U*/A ^ m(l + 2(3/a). (7) 

where n, and m and are integers, while k is a large real 
number, proportional to L, whose value depends on de- 
tails of the system. Since Nq is very sensitive to details, 
it may be considered a random number, different for ev- 
ery nanotube. If a/ (3 and TVo satisfy the conditions (i)- 
(iii) we obtain an infinite sequence of PS electron entries. 
Conditions (i) and (ii) may be slightly relaxed when a fi- 
nite SU is included; in order to find a finite PS sequence 
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in that case the actual U c in Eq. (|l|) has to be in a vicinity 
SU of U* given in Eq. (0), and condition (i) needs to be 
satisfied with an accuracy of order SU/ A. However, finite 
SU makes the condition (iii) somewhat more difficult to 
fulfill. Even with the inclusion of a finite J it is possible 
to find, as shown in Fig. ||g, a set of parameters that give 
a PS sequence. 
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FIG. 2. Evolution of the total spin S, under a non-uniform 
potential, whose maximum value on the nanotube is cV g . Var- 
ious cases are shown, with SU = 0.2A and with other param- 
eters as given in the inset. The total spin S = 0, 1/2 or 1, is 
shown by the thin solid lines. Tick marks on the thick hor- 
izontal lines show the entry points of electrons. The a>axis 
is the product /3cV g , in units of the level spacing A. A small 
magnetic field h — 0.001A was included. 

To illustrate how the PS scenario works, we have cal- 
culated numerically the addition spectrum and the total 
spin of the ground states as a function of the maximum 
voltage on the nanotube, cV g , that varies with the poten- 
tial on the gate electrode. Results are shown in Fig. |^ for 
few selected sets of parameters, which are listed in the 
insert. Case (a) corresponds to an ideal situation when 
all electrons enter with the same (up) spin; in case (c) 
all the electrons enter with spin down. In case (b) No is 
in-between the ideal matching conditions and the spin se- 
quence is alternating. These three examples have J = 
and parameters corresponding to n = in condition (@)ii. 
Cases (d)-(f) correspond to n — 1. We see that by chang- 
ing slightly the model parameters we can get PS series, 
alternating spins, and transitions between them. Finally 
in case (g) we show that even with finite J it is possible 



to get a PS sequence. 

We note that by increasing the value of h by an amount 
of order A, one can convert a sequence of up spin entries, 
as in (a) , to a sequence of down-spin entries as in (c) [|l2| . 
This is consistent with experimental results in I. 

In Figure 1, of I, internal transitions of the nanotube 
appear as discontinuities of slope in the boundaries of 
the regions in the plane of bias voltage and gate voltage 
where conductance is inhibited by a Coulomb barrier. In 
our model, the relative change in slope on transition from 
an A-level to a B-level should be a/ (3, or about 7%. In I, 
however, the changes observed are much larger than this. 
Such large changes may be more easily explained if there 
are localized states at the Fermi energy, due to disorder, 
which could then respond quite differently to the gate or 
bias voltages [n3|. On the other hand, slope changes ob- 
served in Ref. jSj], for a different nanotube, were smaller, 
of order 20%, and might possibly be consistent with our 
model of a perfect nanotube. 

In conclusion, we presented a model of the armchair 
carbon nanotube in a nonuniform external potential and 
elucidated a microscopic mechanism for internal transi- 
tions in which spin flips occur. This mechanism could 
lead to electrons entering in parallel spin sequences, as 
was observed in I, and is consistent with the lack of these 
in II. We note that a sequence of parallel spin entries 
occur when few special conditions are fulfilled, in most 
cases one or more are not fulfilled and there is alternate or 
irregular spin sequence. However, there is a finite prob- 
ability, that could be increased with the proper choice of 
gate configurations (that control the model parameters), 
to find a long sequence of parallel spins entries. 

We highly acknowledge discussions with A. Bezryadin, 
M. Bockrath, C. Dekker, C. M. Lieber, H. Park and 
T. H. Oosterkamp. This work was supported in part by 
the NSF through the Harvard MRSEC (grant DMR 98- 
09363), and by grants DMR 94-16910, DMR 96-30064, 
DMR 97-14725. KB was supported by the Foundation 
for Polish Science (FNP). 

* On leave from: Institute of Theoretical Physics, Warsaw 
University, ul. Hoza 69, 00-681 Warszawa, Poland. 
[1] S. Iijima, Nature 354, 56 (1991). 

[2] N. Hamada et al., Phys. Rev. Lett. 68, 1579 (1992); 

R. Saito et al., Appl. Phys. Lett. 60, 2204 (1992); 

J. W. Mintmire et al, Phys. Rev. Lett. 68, 631 (1992). 
[3] M. Bockrath et al, Nature 397, 598 (1999). 
[4] M. Bockrath et al, Science 275, 1922 (1997). 
[5] S. J. Tans et al, Nature 386, 474 (1997). 
[6] A. Bezryadin, A. R. M. Verschueren, S. J. Tans, and C. 

Dekker, Phys. Rev. Lett. 80, 4036 (1998). 
[7] S. J. Tans, M. H. Devoret, R. J. A. Groeneveld, and C. 

Dekker, Nature 394, 761 (1998), experiment I in the text. 
[8] D. H. Cobden et al, Phys. Rev. Lett. 81, 681 (1998), 

experiment II in the text. 
[9] H. W. Postma, Z. Yao, and C. Dekker, preprint, LT22 

Session llaC. 



4 



[10] Entering of consecutive electrons with parallel spins could 
be explained if the ground state spin of the nanotube was 
large, but it can be ruled out because (i) bulk graphite 
layers are not spin polarized and (ii) it does not give an 
explanation for the different in the experimental results 
of I and II. 

[11] P. Fulde, Electron Correlations in Molecules and Solids 

( Springer- Verlag, Berlin - N.York, 1995). 
[12] K. Byczuk, Y. Oreg, and B. I. Halperin, manuscript in 

preparation. 

[13] N. S. Wingreen and A. Odintso v, March meeting: BAPC 
[JC02.08], 1999, L. Balents, |cond-mat/9906032| , 1999 
Moriond Les Arcs Conference Proceedings. 



5 



